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Abstract 

For a general class of saddle point problems sharp estimates for Babuska's inf-sup stabil- 
ity constants are derived in terms of the constants in Brezzi's theory. In the finite-dimensional 
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1 Introduction 

In this paper we consider linear problems in saddle point form: Find u G V and p G Q such that 





u 




~f 




p. 




_8_ 



with 



A B* 
B -C 



(1) 



with complex Hilbert spaces V and Q, f eV* and g £ Q* and linear operators 

AeL(V,V*), B G L(V,Q*), B*eL(Q,V*), C G L(Q,Q* 
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Here H* denotes the dual of a Hilbert space H, which is defined as the set of all bounded and 
antilinear functionals on H, L(H\,H2) denotes the set of all bounded and linear operators from 
H\ to H2, and B* is the adjoint of B. 

For the special finite dimensional setting V = C" and Q = C m , the linear operators above can 
be identified with corresponding matrices 



BeC" 



B* e C 



CeC 



where B* denotes the conjugate transpose of a matrix B. Then (OQ) becomes a linear system of 
equations in standard matrix- vector notation. 

A variety of interesting problems in science and engineering can be framed in this way, for 
example, mixed formulations of elliptic boundary value problems, the Stokes problem in fluid 
mechanics and the optimality system of PDE-constraint optimization problems, see, e.g., lfT8l . 
ffHl ^ El , (261, 021. Nonlinear versions of such problems are typically approximated by a se- 
quence of systems of the form ® during a linearization process, see, e.g, |26| . Discretized ver- 
sions of these problems lead to finite dimensional linear systems of the form (OQ) with large scale 
and sparse matrices A, B and C, see, e.g, (8]|. The use of a framework based on complex rather 
than real Hilbert spaces is motivated by particular applications in spaces of time-periodic func- 
tions, for which the Fourier transform naturally introduces a representation in complex Hilbert 
spaces. 

We will address two fundamental and related topics: stability estimates for (OQ) in general 
Hilbert spaces and spectral properties of the matrix in the finite dimensional setting for Her- 
mitian matrices ^ . 

Stability estimates for (OQ) are estimates of the norm of the solution 



ex=VxQ 



in terms of the norm of the right-hand side 

\f 



G V* x Q*. 



The product space V* x Q* can be identified with X*, the dual of X, in a canonical way. That 
makes ^ a bounded linear operator from X to X* and a stability estimate can be written in form 
of an upper bound for the corresponding norm of ^# : 



-1 1 



\L(X*,X] 



< 



(2) 



where || • denotes the norm in a Hilbert space H. Of course, such an estimate depends on the 
choice of the norm or better the inner product in X. We will concentrate on inner products in X 
of the particular form 



(x,y) x = (u,v)v + (p,q)Q for x ■■ 
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V 


X = 




, y = 






p 


9 



(3) 
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where (•,•)# denotes the inner product in a Hilbert space H. 

Stability estimates are required for showing the well-posedness of (QQ) and are also important 
for discretization error estimates. Rather sharp stability bounds c for the case C = are already 
contained in the pioneering paper [7] and have been improved, e.g., in [|27l . For C = we will 
show how to further improve the estimates and present the best possible bounds in a general 
Hilbert space setting under the same assumptions as in [7] and Il27ll . For general linear operators 
C we will formulate a simple but rather helpful criterion which leads to stability estimates for 
certain problems from optimal control. 

Note that, contrary to ©, sharp upper bounds of the form 

\\^L(x,x*)<c (4) 

are easy to obtain and well-known. 

For V = C" and Q = C"\ inner products can be represented by Hermitian and positive definite 
matrices. In particular, © corresponds to inner products of the following form 

(x,y) x = (Pu,v) + (Rp,q) = (^x,y) with & = 

where (-, •) denotes the Euclidean inner product in C and P G C nxn , R e C mxm , and, therefore, 
also G £,( n+m ) x ("+'") are Hermitian and positive definite matrices. 

In the Hermitian case = the estimates © and © immediately lead to corresponding 
estimates 

c < |jU| <c, 

for the eigenvalues \x of the matrix M = \M ' . Therefore, we have 

A* e [-c,-c]u[c,c]. 

We will show that inclusions of this form with two intervals symmetrically arranged around 
are appropriate for the case C = A and R = P due to symmetry properties of the spectrum of <Jt . 
For the case C = and A positive definite, sharp estimates of the more general form 

jl e [/il,jU2]U[/i 3 ,jU4] (5) 

with Hi < H2 < < /i3 < H4 were derived in E4ll . For positive semidefinite C and indefinite A, 
rather sharp estimates were presented in lfT2l . For the special case C = we will further improve 
these estimates and present the best possible bounds under similar assumptions as in lfT2l . 

Preconditioned Krylov subspace methods are an important class of iterative methods for solv- 
ing (OQ) in the finite dimensional setting with large and sparse matrices A, B, and C; we refer to 
the survey article [4] and the references cited there for these as well as other solution techniques 
for saddle point problems. In the Hermitian case ^* = ^ a prominent representative of these 
methods is the preconditioned minimal residual method (MINRES), see 11231 . whose convergence 
analysis relies on inclusions of the form ©, see, e.g., E3l . |[T3~1 . 

The paper is organized as follows. Sect. |2] contains stability estimates of the form © in 
general Hilbert spaces for C = 0. In Sect. [3] we study the finite dimensional Hermitian case and 



P 
R 
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derive sharp estimates of the form © for C = 0. Structural properties of the spectrum of for 
C = A and R = P are studied in Sect. @J In Sect. [5]the general results of the preceding sections are 
applied to two model problems from optimal control, distributed time-periodic parabolic control 
and distributed time-periodic Stokes control. The paper finishes with concluding remarks and 
some technical details of a proof shifted to the appendix. 



2 The general case 

For the analysis we use a reformulation of £[]) as a variational problem: The linear operators A, B, 
and C uniquely determine sesquilinear forms a, b, and c on V x V, V x Q and Qx Q, respectively, 
given by 

a(u, v) = (Au, v) , b(y, q) = (Bv, q) = (B*q, v) , c(p, q) = (Cp, q) , (6) 

where, in general, (•,•) denotes the duality pairing in Hilbert spaces. (It reduces to the Euclidean 
inner product, as used in the introduction, in matrix-vector notation in the finite dimensional 
setting.) Then (Q~|) becomes a mixed variational problem: For given / 6 V* and g E Q*, find 
u EV and p e Q such that 

a(u,v)+b{v,p)=f(v) forallvGV, 
b(u,q)-c(p,q) = g(q) for all q G 2, 

or, equivalently, a variational problem in the product space X = V x Q: Find x EX such that 

3§(x,y) = &(y) for all y e X (8) 

with 



(x,y) = (^x,y)=a(u,v)+b(v,p)+b(u,q)-c(p,q), &(y) =f(v)+g(q) 



for 



x - 



and y 



From the famous Babuska-Brezzi theory, see O, J3), 0, applied to the case C = 0, it is 
well-known that the following three statements are equivalent: 



1 . .J^f is an isomorphism between X and X* 

2. < oo with 



K*,y)\ 



sup sup 

O^xeXO^yeX \\X\\X\\y\\X 



and 



inf sup 



\&(x,y)\ 



inf sup 



o^xexo^yex \\x\\x\\y\\x o^yexo^xex \\ x \\x\\y\\x 



= r>o. 
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3. ||a|| < oo, \\b\\ < °° with 

sup sup 

( 

and 



\a(u,v) 



inf 



O^ueVQ^veV ll M lln|v||y 
a(u,v)\ 



sup 



inf 



\b(v,q)\ 

SUp SUp Tj Tj — .. .. 

\a(u,v)\ 



sup 



a > 0, 



l*(v,9)l 



inf sup 



j8 >0. 



It is also well-known that the quantities 
associated operator 



o^qeQOytvev \\ v \\v\\q\\Q 

m and 7 can be expressed in terms of norms of the 



1 

r 



(9) 



\L(X,X*)i ~ - ||^~ I If A' .Xr 

With these notations and relations, © and (HJ) take the following form 

0<c<7 and \\3§\\<c<°°. 

Our aim is thus to derive lower bounds of y and upper bounds for . One possible approach 
to obtain these bounds is based on available lower bounds for a, /3 and upper bounds for ||a||, 
\\b\\. These quantities can also be expressed in terms of norms of the associated operators A and 
B. We follow here the presentation in [1J and introduce a space decomposition on V: 

V = V + Vi with V = kerfl and V\ = Vq", 
and define the operators Aij G L(Vj, V*) for all /, j G {0, 1} by 

(AjyK,v) = a(w,v) for all w g V/, v G Vi- 
and Bi €L(Vi,Q*)by 

(fiiv, <?) = b(v, q) for all veVi, qeQ. 
Then it is easy to see that Aqq and Bi are isomorphisms and we have: 



L(V,V*), 



\B\ 



L(V,Q*) 



\B 



and 



1 



1 



- - \\ a qo\\l(v*,v )i jg-H 5 ! \\L(Gr,v l ) = \\V 3 i) llL07.fi)- 



HIl(Vi,2*) - 
5 *w 



l 5 lllL(fi,^ 



(10) 
(11) 



With these operators we can rewrite the original problem in a 3-by-3 block form, as suggested 
in El: 



Aqq Aoi 
Aio A u 5* 
Bi 



"0 




fo 


Ml 




fl 






g 



(12) 
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with u = uq + u\, Uj G Vi and// G V,* are given by (fi,v) = (/, v), v € V - / for / G {0, 1}. As observed 
in [OQ the inverse of the 3-by-3 block operator in (PT2l) is given by 



r 1 

1 oo 




AiqA 



-l 

00 








-A OQ A iB 1 1 



(S^'^il-AioAoo'Ao^B,- 1 



Using this representation and (fTOl) . (fTTT) it follows immediately for (fT2l) that 



j8 II^OO^ 01 ! 



"IKIk" 




r i o 


ll»illv 


< 





.bile. 







jff ^ || A n _A io A oo A oi|| 





'Wfohf 




ll/illv? 




.klle*. 



Here, we dropped indices of the operator norms for simplicity. It is clear from the context which 
operator norm is meant. Obviously, we have ||A ;; || < ||A|| = ||a|| for i G {0, 1}. Therefore, 



A oo A 01 II < 



^00 



a 



AiqA 



-l I 

00 I 



< 



and 



which imply 





An- 






1 "ll W 


< 


.bile. 





An -AioA^AoiH < ||a|| + 



a 

|2 
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a 



f 



HI 1 HI [ii HI 
a/3 p ^2 V « 



0\\V* 



1 1 1 V,* 



.klle* 



Using 



IM|v < HkoIIv + ||«i||v an( * 
it easily follows that 



0W 



< 



IV," 



< 



||m||v 


< 


.bile. 





1 

a 

1 (a , IN 



1 M _i_ H 
l i+ a 



/3 2 



llslle* 



(13) 



(14) 



(15) 



These are exactly the estimates that can be already found in and later in [Q]| with essentially 
the same line of arguments. Since 
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~f 


p. 




g 



and / G V*, g G Q* can be chosen arbitrarily, this estimate shows that 
/_1 Hl(x*,x) <P(Di) with D x - 



1 ( i , « 



I n H \ f 1 _i_ H 

j6 \ + a ) p l 1+ a 



(16) 
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where p (M) denotes the spectral radius of a matrix M. From ® we learn that 1 / p (Di ) is a lower 
bound for 7. 

With a slight variation of the arguments we obtain a better lower bound for 7. The essential 
step is an improvement of the estimates in (fT3l) . (fT4l) . 



Lemma 1. Wfe /save 



IIA10A 



00 



< 



a 



, -1, ||A^Aoi||< 



-1 



a- 



lAn -AioAJAqiII < 



Proof. By using the simple estimate 



it follows that: 



I + IIA10A 



A o 
Mo 



-1m2 
00 II 



< 



00- 



Aoo ^01 

A10 An 



A10A 



-1 
00 . 



Aoo 
Mo 



r 1 

^00 



< 



a 



2 ' 



from which the first inequality follows immediately. The second inequality can be shown analo- 
gously. For the third estimate we have 



IIAioAjAoi - An 



[AioAoo 1 -/] 



A01" 




An. 





<||[AioA^ -/] 



A i 
An 



l + l|AioA 00 J 



1||2 



A i 
An 



\\a\\ 

< \\a\ 

a 



a 



which completes the proof. 

Using the estimates of the previous lemma we obtain the following result. 
Theorem 1. With the notation introduced above, the following holds: 

1. Let y opt (a,P, ||a||) be the smallest positive root of the cubic equation 

Ai 3 -(IH| 2 + /3 2 )Ai + a/3 2 = 0. 

Then 

2. We have 

\\u\\v 

< 



□ 



(17) 



|M||y 

\p\\q 



1 

a 



iM 

j8 a 

iM H H 

/3 a /3 2 a 



V* 

g\\ar 



(18) 



and y> a/(l + k 2 ) with k = ||a||/jS. 
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Proof. From (PT51) and Lemma Q] it follows that 



"||"0 


\v 




ll/o 


\v*~ 






a 





| Ml 


\v 


<E 


Wfi 


\v* 


, with E = 










JH 


Q. 




.\\g\ 


0* . 








1 

J8 



1 / a 
1 

jS 2 a 



-1 



(19) 



Therefore, 



<P(E) 



x 



f 



X* 



which implies that l/p(E) is a lower bound for y. Since E is a nonnegative matrix, p(E) is an 
eigenvalue of E, see E page 26]. Therefore, p (£) is the largest positive root of the characteristic 
polynomial of E, which directly leads to the cubic equation (fTTT) for \/p(E) = 7 op t(«,/3, ||a||). 
Using 



,|2 || i|2 
\ u o\\v + ll w l II V i 



^ = nybiife +11/1 life 



and Cauchy's inequality, it easily follows from (1191) that 



|| M || V 




"ll/llv*~ 






<D 2 




, with D2 


He 




llsllfi* 





1 in 

IM H J!fL 

/3 a j8 2 a 



As above, this implies that 1 / p (D2) is a lower bound for 7. Observe that D% is a rank-one matrix 



D 2 = - 
a 



[l k] , with K 



P ' 



whose spectrum is given by {0, ( 1 + K" 2 ) / a}. Therefore, 1/p (D2) = Ct/( 1 + K" 2 ), which com- 
pletes the proof. □ 

We would like to remark that the estimate 7 > 73(05, /3, ||a||) is sharp. Indeed, for the matrix 
Jt in (0Q) with 



a 



-vTI 



VNI 2 -« 2 ' 



all 2 - a 2 



we obtain 



a 



1 

a 





B = [0 j8] and C = 0, 



(20) 



* 



1 * 



or 
/3 

1 M 

/3 1 a 



The conditions of the Babuska-Brezzi theory are satisfied with constants ||a||, a, and /3 for ^ 
/, the identity matrix. Since coincides with the matrix E, it immediately follows that 



7=7o P t(a,j8,||a||). 
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We would also like to emphasize that both bounds of TheoremQ]are sharper than the classical 
bound 1/ p(D\). Indeed, if we apply the estimates for 7 from Theorem \T\ to the matrices in (1201) . 
we obtain 7 = 7 op t(ce, ||a||, \\b\\) > cc/(l + k 2 ). Actually, we have a strict inequality 

cc 

y opt (a,\\a\\,\\b\\) > — -j, 
since the cubic polynomial in (fTTT) is strictly positive on [0, <x/(l + fc 2 )]. Moreover, for 

we have 



1 

a 



1 k£ 

Kt, K 2 ^ 



D 2 =D{\) and D { = D(^) with = "jy. > L 

It is easy to check that p(D(^)) is a strictly increasing function of <§. Therefore, p(D\) > p(£>2)> 
or, equivalently, 

a 1 

r+^pw 

Comparing the results in Theorem Q] with similar bounds in the literature, we notice that the 
bound presented in ll27l is of the same form, say f((X, j8, ||fl||), and one can show that 

a > f(cc, J3,||fl||)> 1 



l + K 2 ' '^" f " 11 — ' p(Di) 
For completeness we finally include a trivial upper bound for 



Since we have 



~\\f\\v*~ 


< 




ii*ir 




||«||v 


JHIe*. 











.bile. 



it follows that 



<P 



|a|| \\b\ 
\b\\ 



« + \ \\a +4p> 



(21) 



It is well known that this estimate is sharp; see, e.g., [24J. 

3 The finite dimensional Hermitian case 

In this section we consider the finite dimensional case V = C n and Q = C m . Inner products 
in these spaces can be represented by Hermitian and positive definite matrices P G C" x " and 

g (^ynxm. 

(u,v)v = (u,v)p, (p,q)Q= {p,q)R, 

where here and in the sequel the following standard notations are used: For a Hermitian and 
positive definite matrix M E C xr , the associated inner product is given by (z, w)m= (Afz, w) for 
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Z, w G C. Both the vector norm and the matrix norm associated with the inner product (•, -)m are 
denoted by \\-\\m- 

With the notation introduced above, the inner product ©inX = Vx<2 = C n+m is given by 



(x,y)x = (x,y)g* with & 



P 
R 



+m) x (n+m) 



(22) 



The discussion in this section is restricted to the Hermitian case: 



For the inf-sup constant y and the norm || 
introduced in Sect. [21 it is well-known that 

7=lMmin| and 



of the bilinear form S3 associated with j$ as 



|M 



max 1 5 



where jj,^ and are the eigenvalues of the matrix = with minimal and maximal 

modulus, respectively. Or, equivalently, for the eigenvalues /! of the generalized eigenvalue 
problem 

Jix = \l ,<?>x, (23) 

we have 

^ [-11^11, -yj u [y,p3||]. (24) 

For the rest of this section we concentrate on matrices of the form 



A B* 
B 



G (Q("+ m ) x 



(25) 



with A* =Ae C nxn satisfying 



(Av, v) > for all ^ v G kerfl, 

and B G C mxn of full rank m < n. Under these conditions the matrix «/# is Hermitian and 
non-singular. This setting is typical for optimality systems (Karush-Kuhn-Tucker conditions) 
of equality-constrained optimization problems. Let Pi G L(V,-, V*) be given by 

(Pjii, v) = (Pu, v) for all u, v G V;, i G {0, 1}. 

Under the conditions above we have the following well-known alternative representation of the 
inf-sup constant a: 

(A 00 vo,v ) 



a 



inf 



(26) 



O^voeVb (PoVo,vo) ' 

Instead of ||a||, the norm of the bilinear form a associated with A in V, we assume more detailed 
information on A in terms of the extreme eigenvalues of P l A, which can be written as 



lA 

Amin o^vev(Pv,v) 



Av,v 

SU P 75 — v 

O^vev (Pv,v) 



(27) 
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Of course, we have A^ ax > a > and ||a|| = max{|AAj, A max }. 

The following representations for the inf-sup constant /3 and the norm \\b\\ of the bilinear 
form b associated with B in Q hold: 

2 . (BP l B*q,q) . ^m, vQ 

p = inf — : = inf — ^ -— : , (28) 

o^qeQ {Rq,q) Q^v x eV x (PiVi,Vi) 

and 

„ (BP 1 B*q,q) 

\\b\\ 2 = sup K - fit- . (29) 

Since ^# is Hermitian and indefinite we have an inclusion of the form 

with jUi < jU2 < < jl?, < /I4 for the eigenvalues /i of the generalized eigenvalue problem (1231) . 

Sharp bounds jUi, jU2, and ^4 in the case AA n > can be found in ll24ll . In lfT2l these bounds 
were extended to the case AA < and read, in general, 



1 



Ml = 2(Cn-V(Cn) 2 ^IHP 
1 



M2=2(*max-\/(^ m ax) 2 +4/3 2 ), (30) 



1 



M4=2^max+V( A max) 2 +4||^|| 2 

It remains to discuss lower bounds for the positive eigenvalues. For the case AA n > a simple 
bound is known (see E4ll ): 

M3 = ^min- 

This bound is sharp for problems with a = AA . Rather sharp bounds with no restriction on A A 
can be found in [fT2l . The best bound presented in lfT2l is the smallest root of a cubic equation 
whose coefficients are given in terms of a, AA n , A^ ax , and J3. For the case AA n < we will 
derive now a sharp bound that involves a similar cubic equation. 



Theorem 2. Assume that AA n < and let jl be any positive eigenvalue of 8? \M. 
1. Let y pt((X,P i^mini ^max) ^ e ^ w smallest positive root of the cubic equation 

M 3 - {*L + AiJ M 2 + {*L*L - /3 2 ) M + «/3 2 = 0. (31) 
Then we have jU > y opt (a, j8, AA n , A A J 
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2. We have 



ju > < 



X A . X A + B 2 
/3 2 



3 A 3 A 

min max 



2(X A . +X A ) + \ V 2(A A . +X A 

V min 1 max) \ \ V^min 1 mew 



mm max r 1 



+ 



l A +l A 

min 1 max 



Proof. Let /i > be an eigenvalue of the eigenvalue problem (1231) with eigenvector 
x = [uq, ui,p] T ^ 0. Using the framework introduced in [0Q, see Sect.[2l this eigenvalue problem 
reads 



A o 


Am 


0" 




«o 




>o 





0" 




wo 




An 


B\ 




mi 


= M 





Pi 







Ml 





Si 

















7? 







i.e., 

Aoomo+Aoimi =nP u , 
AioMo+Anwi +B\p = iiP\u\, 
B\u\ = \xRp. 

First we consider the case that Aoo — /i Pq is non-singular. Then we obtain from the first and the 
third equations 

m = -(Aoo-juifo) _1 Aoi"i and p = —R 1 B i uu 

which immediately implies that u\ ^ 0. Using these relations for eliminating uq and p from the 
second equation and taking the inner product with jiiui we obtain 



fi 2 (Piwi, wi) +ju (Aio(Aoo - M A)) ^oiwi^i) 

-/l (AiiKi,Ki) - (5ii? _1 5ii<i,i/i) = 0. 



(32) 



We have (Aoovo, vo) > oc (Pqvo,vo) for all vo € Vb; see (|26l) . Then, for < /i < a, it is easy to 
check that 



((Aoo-fiPo) Vvo) < 



a 



(Aqq 1 v , v ) for all v € V • 



a-ji 



With vo = Aoi«i, an upper bound for the second term on the left-hand side in (1321) follows: 
H (A w (Aoo - nP y l A 0l ui,ui} < (AioA^)Aoitti,«i)- 

0£ jLi 



For the last term on the left-hand side in (1321) we have 

{B\R- l Biui,ui) > /3 2 (Pimi,mi), 



(33) 



(34) 
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see (Eg). Using ([34]) and ([33]) we obtain from ([32]): 



/I 2 (PlU U Ui) + 



/MX 



(aJ.Aqiu\,Aq\u\) -/i (An ui, «i) -/3 2 (Piwi,mi) > 0. 



After dividing by (Piwi, «i) it follows that 



a 



a — ji 



r\-r 2 



/3 2 >0, 



(35) 



with the Rayleigh quotients 



{A\qA^Aq\U]_,u\) 



and r2 



(Aii«i,mi) 



(Pi«i,«i) (PiUum) 
One can show that (see the appendix for the technical details): 



a-/i 

Then it follows from (1351) that 



a ^ (^max + ^min ~ CL )\ L ~ ^max^ 



A 

min 



M 2 + M 



a — n 



(^max ^min ^0 M ^max^n 

a - jU 



(36) 



-j8 2 > 0, 



i.e. 



= M 3 - (Ai x + XLn) M 2 + (Ai x A^ in - /3 2 ) M + a/3 2 < 0. 



Therefore, /i cannot lie in the interval between and the first positive root, denoted by 
7o pt (a, /3, A^ in , A^ ax ), of the cubic polynomial q(fi), because there the cubic polynomial is strictly 
positive. 

If Aoo — fiP is singular, then /i > a because of (1261) . Since q(tt) < 0, it follows also in this 
case that /i > 7 opt (a,/3, A^ in , A^ ax ). That completes the proof of the first part. 

The proof of the second part follows the same line of arguments as presented in the proof of 
Proposition 2.2 in lfT2l and is omitted here. □ 



We remark that the estimate /i > y op t(a,/3, A^ in , A^ ax ) is sharp: For the matrix M in (125 



with 



a 



(*4x-a)(a-*tain) 

A , 1A 



1 A i 1 A p, 
A max ' ^min 01 



'(^ ax -a)(a-A^ n ) 
and = /, we obtain the following characteristic polynomial 

M 3 - (^-max + ^min) + (^taax^min ~ fi 2 ) H + a P 2 , 



B=[0 P] 
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which coincides with the cubic polynomial in (fTTT) . Therefore, the smallest positive eigenvalue 
of J( is equal to y opt (a,/3, A^ in , A max ). 

In many applications A is positive semidefinite with a non-trivial kernel. Then A min = and 
^-max = ll a ll' which leads to the cubic equation 

H 3 - \\a\\ n 2 - p 2 n + a p 2 = 0, 

whose smallest positive root is a sharp lower bound for the positive eigenvalues \l of \M 
in this case, see Theorem [2] and the example above. The second part of Theorem |2] yields the 
following simpler bound: 



2 IMI V V / /3 + v / /3 2 + 4a||a|| 

for the positive eigenvalues }i of 2?~ x ^t '. 

Finally, we would like to stress that because of simple monotonicity arguments, all presented 
estimates on the spectrum of £?~ x ^, remain valid if a, j8, y, and A min are replaced by lower 
bounds and A^ ax , \\b\\, and \\^\\ are replaced by upper bounds. For example, if A^ in is replaced 
by its lower bound — ||a| and A max by its upper bound \\a\\ in (1311) . then we also obtain a lower 
bound for the positive eigenvalues of which is not necessarily sharp but it goes without 

specific knowledge on the spectrum of A other than the spectral radius of P A. Observe that, 
with these replacements, (I3TT) coincides with (flTT) . and the corresponding smallest positive root 
is identical to y op t(a,/3, ||a||), derived in Sect. 12 



4 On a class of matrices with symmetric spectrum 

In this section we specialize our considerations to the following Hermitian matrix 



A B* 
B -A 



eC 



2nx2n 



with A E W xn real and symmetric positive definite, and B e C" x " complex symmetric, i.e., 
B = B T , where B T denotes the transpose of the possibly complex matrix B. We know that the 
matrix ^ has all real eigenvalues, n positive and n negative ones. We next show that the negative 
eigenvalues are the mirrored images of the positive eigenvalues. This property has a few conse- 
quences, both in the choice of the preconditioner, and on the convergence of MINRES. Indeed, 
the spectrum of ^# is symmetric with respect to the origin, and MINRES behaves like CG on a 
matrix having only the positive eigenvalues, but with twice as many iterations. Therefore, MIN- 
RES on will only make some progress every other iteration, showing complete stagnation 
otherwise; we refer to [fTTTl for a similar phenomenon for 2 x 2 block matrices with a different 
nonzero structure. 

We first need the following technical lemma. 
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Lemma 2. Let H be a nonsingular complex symmetric matrix ( i.e., H = H T ), and S be a complex 
skew -symmetric matrix (i.e., S = —S T ), both of size n. Then the (nonzero) eigenvalues of the 
In x In matrix 

"0 I 
H S 

come in pairs, (jU, — jli). 

Proof. For nonsingular H, we have the similarity transformation 







"0 / 




m r -2 




-iH^ 


I 




H S 




/ 




ml s 



where the matrix in the right-hand side is complex skew- symmetric. Since any skew-symmetric 
matrix T is similar to its transpose T T = —T, to each Jordan block with eigenvalue /! in the 
Jordan decomposition of T there is a corresponding Jordan block of the same size with eigenvalue 
-ll. □ 

Theorem 3. Assume thatB is nonsingular. Then the eigenvalues /i of \M come in pairs, ([1, — [l), 
with \x Gl. 

Proof. The eigenvalue problem for ^# can be written as 

Ax + B*y = /lx, Bx-Ay = }ly. 
Substituting x = B ' ([ll +A)y in the first equation and after some simple algebra we obtain 

-ll 2 B- l y + /i [AB 1 - B l A]y + [AB l A+B*]y = 0. 
After multiplication by B^ from both matrix sides we obtain 



-\L 2 y 



B^AB 3 -B 3A53 



y + 



B l iAB l AB l i+ BW*B? 



with y = B iy, which gives 

-ll 2 I + /l(G-G T ) + (GG T +B?B*B?) 



y = 0. 



The eigenvalues of the quadratic matrix equation above can be obtained as eigenvalues of the 
linearized problem 



O I 

GG T +B?B*B? G— G T 



Z = jUZ & &Z = jUZ. 



We then recall that all eigenvalues of and thus of , are real, and that n of them are positive 
and the other n are negative. The matrix ^ satisfies the hypotheses of Lemma [2l therefore the 
eigenvalues of come in pairs (jU, — /x), which completes the proof. □ 
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5 Application to optimal control problems 

Let Q. be an open and bounded domain in W 1 for J e {1,2,3} with Lipschitz-continuous bound- 
ary r. For T > 0, we introduce the space-time cylinder Q T = £1 x (0, T) and its lateral surface 

z r = rx(o,r). 



5.1 Distributed optimal control for time-periodic parabolic equations 

First we consider the following model problem: Find the state y(x,t) and the control u(x,t) that 
minimizes the cost functional 

J{y^ u ) = \ I I \y{x,t)-y d (x,t)\ 2 dxdt + ^r f f \u(xj)\ 2 dxdt 
2 JO JO. 2 Jo Jo. 

subject to the time-periodic parabolic problem 

d < 



Ay(x,t) 


= u(x,t) 


in Q T . 


y(x,t) 


= 


on 


y(x,0) 


= y(x,T) 


on 


u(x,0) 


= u(x,T) 


on £1. 



Here yd(x,t) is a given target (or desired) state and v > is a cost or regularization parameter. 
We assume that yd(x, t) is time-harmonic, i.e.: 

2nk 

yd(x,t) = yd(x) e l03t with (0 = —^- for some keZ. 
Then there is a time-periodic solution to the original control problem of the form 

y(x,t) =y(x)e iat , u(x,t) = u{x)e iat , 
where y(x) and u(x) solve the following time-independent optimal control problem: Minimize 

\ I \y{x)-y d {x)\ 2 dx+^r I \u{x)\ 2 dx 
2Jei 2 Ja 

subject to 

i(Oy{x) — Ay(x) = u(x) in 

y(x) =0 on r. 

Using an appropriate finite element space V/ 2 of dimension n for both y and u, we obtain the 
following discrete version: Minimize 

\(y-yJ*M(y_-y_ d ) + ^u*Mu 
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subject to 



i(0My + Ky =Mu. 

Here the real matrices M and K are the mass matrix, representing the L 2 -inner product in Vh, and 
the discretized negative Laplacian, respectively. The underlined quantities y, y,, and u denote 
the coefficient vectors of the corresponding finite element functions relative to the chosen set of 
basis functions in Vh- 

The Lagrangian functional for this constrained optimization problem is given by 

■&(y, u,P) = \ (l -hY M (l ~ l d ) + \ M-* M <± + f O m Z +Ky- Mu) , 

where p denotes the Lagrangian multiplier associated with the constraint. The first-order nec- 
essary optimality conditions, which are also sufficient for the problem considered here, are 
VJ?(y,u,p) = 0, and read in details: 

M K-icoM 
vM -M 
K + icoM -M 

This is a linear system of the form (OQ) with 

eR 2nx2n , B=[K + icoM -M]eC nxZn , C = 0. (38) 

The system (1371) was discussed in G5Tl for the special case £0 = 0, which corresponds to an elliptic 
optimal control problem. Observe that in this case all matrices are real. The preconditioner 
constructed in ll25l is an indefinite 3-by-3 block matrix and leads to convergence rates of the 
preconditioned conjugate gradient method which do not deteriorate if the mesh size, say h, and/or 
the cost parameter v approach 0. 

Based on ideas very close to those developed in E51l we obtain the following properties of 
the bilinear forms a and b associated with A and B, respectively, in terms of the quantities a, 
^min' ^max^^ and \\b\\: 

Theorem 4. Let ^# be given by (0) and ( TJffl) . Then, for given by 

, R = ^Y, andY =M+^v(K+coM) 1 
we have 

a>2-V2, A,4„>0, Ai x <l, j8 > ^, \\b\\ < 1. 

Proof. The estimate A^ in > is trivial. The upper bounds for A^ ax = \\a\\ and \\b\\ follow com- 
pletely analoguously to the proof of Lemma 4.1 in ll25l . 



y 




My, 


u 







p 








(37) 



M 
vM 



P 
R 



with P = 



Y 
vM 
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The proof of the lower bounds for the inf-sup constants a and /3 in ll25Tl covers only the case 
CO = 0. An essential step of that proof was the estimate 



K 



< 



\m \\u\\m 



for all y and u satisfying the state equation for (0 = 0, i.e., Ky = Mu. That part of the proof has 
to be replaced for general CO by the estimate 



\l + co\ 



\m < v^HvIImIIwIIm 



for all y and u satisfying the state equation (K + icoM)y = Mu, which easily follows: 

(llj|||+«lblllf) 2 <2(lly|||+« 2 llj|ia 



2 \ ((K + ia>M)y,y)\ 2 = 2 \(Mu,y)\ 2 < 2| 



\M\ 



\M- 



All other arguments are completely identical to the corresponding arguments used in Il25ll and 
are omitted. □ 

If a, A^ in , A^ ax , /3, and \\b\\ in (1301) and (1311) are replaced by the corresponding lower or 
upper bounds provided by Theorem |H it immediately follows that the spectrum of ^ is 
contained in the set 

1 



u 



jU3,2(l + V5 



where ji^, is the smallest positive root of the cubic equation 

A i 3 - i u 2 -- i u + l- — = 0. 
These intervals read in 3 -digit accuracy 

[-1, -0.366] U [0.396, 1.618]. 
From the second part of Theorem [2] we know a simple lower bound for ^3: 

4-2V2 



M3 > 



1 + VlV-Sv 7 ^ 



0.346. 



There is a frequently used alternative approach for solving (|37l) . From the second row of 
(1371) it follows that u = p/v. Therefore, we can eliminate the control u and obtain the reduced 
optimality system 

|" M K-icoM] 
K + icoM -i. 

This system was discussed in [28J for the special case co = 0. The preconditioner constructed in 
Il28l leads to convergence rates of the preconditioned MINRES method which do not deteriorate 
if the mesh size h and/or the cost parameter v approach 0. The results from [|28l were extended 





y 




My, 




p 
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to the case CO 7^ in [fTTTl (based on results from fTOlO. where a preconditioner was constructed 
and a bound for the number of MINRES-iterations was derived which is independent of h, v, and 
CO. 

Here we will shed some new light on the preconditioner from ifTTl by presenting a slightly 
different analysis of the already known properties within the framework of complex matrices, 
and by supplementing these properties by a new statement on the symmetry of the spectrum of 
the preconditioned matrix. This new approach is helpful for extending the analysis of precondi- 
tioners for other optimal control problems, like the one discussed in the subsequent subsection. 

A simple scaling leads to the equivalent system 

M y/v(K-iG)M) 
^(K + icoM) -M 

which is of the form (OQ) with 

A=MeR nx '\ B = y/v(K + icoM) eC nx '\ and C = -A. 

Using Sect. [2] and |4] we obtain the following results on the eigenvalues of 

Theorem 5. Let Ji be given by (0) and ( TJ91) . Then, for s , given by 





y 
















(39) 



P 
P 



with P real and symmetric positive definite, the spectrum of 1 jtfi is real and symmetric around 
zero. Moreover, for P = M + y^v (K + COM), the following estimates hold 



7> —j= and 
v3 



< 1, 



where y is the inf-sup constant, and 3$ is the bilinear form associated with ^ . 

Proof. The symmetry of the spectrum of around zero directly follows from Theorem[3] 

applied to the similar matrix ' . Let 

/ (1-/)/ 
(1+0/ -/ 



Then, by direct calculations, one shows that 

Re(^#jc, Jrffx) = \\x\\ 2 gg and \\Jifx 

Therefore, 

. „ \(j?x,w)\ 
7 = inf sup 



2 

\9> 



311x111, foralljceC 2 ". 

\{J(x,^x)\ 



(u.v4.y(Um-€A' l- v l ^ll 11 '!!-^ o^xex \\A\ &>\\^ x \\ 



> inf Re (^ x '^ x ) 



\9> 



inf 

O^xeX ||x|| r ^V3||x||,$ja 



1 

7T 
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Moreover, since 



(j£x 7 x) = (My,y) + 2y/v(Re(Ky,p) - a>\m{My,p)) - (Mp,p) for x 



and 



\Jm(My,p)\ < \(My,p)\ < MmWpWm, \MKy,p)\ < \(Ky,p)\ < IMHHk 
it follows that 



'x,x)\ < \\y\\i + 2Vy(\\y\\K\\p\\K + co\\y\\ M \\p\\M) + \\p\\ 2 M 

< lbllM + V / v(l|j||l+INll + ® (IMlM + lbllAf)) + HP| 

*x,x) for all ;t e C 2 '\ 



2 
M 



which implies that 



x,x) 



\(J?x,w)\ 
sup sup -pj — j-: — - — - — = sup 



< 1. 



This completes the proof. 



□ 



If 7 is replaced by the lower bound l/v3 and by the upper bound 1, it immediately 
follows from d24l) that the spectrum of 3P~ X ^( is contained in the set 



-1. 



1 

V3 



U 



[-1, -0.577] U [0.557,1]. 



Remark 1. The essential step for estimating the inf-sup constant 7 from below was the introduc- 
tion of the matrix M '. Translated into a general Hilbert space setting and slightly more general, 
the essential requirements on such a linear operator Jrf? : X — > X are that there are some positive 
constants c\,c 2 such that 

|^(z,jrz)|>ci||z||i and || < c 2 \\z\\x for all z eX. 

Then it follows analogously to the previous proof that 7 > c\jc%. 

5.2 Distributed optimal control for the time-periodic Stokes equations 

Next we consider the following problem: Find the velocity u(x,?), the pressure p(x,t), and the 
force f (jc, t) that minimizes the cost functional 

7(u,f) = ^/ / \u(x,t) -u d (x,t)\ 2 dx dt+^- [ [ \f(xj)\ 2 dxdt 
2 Jo Jq. 2 Jo J & 
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subject to the time -periodic Stokes problem 



—u(xj)- 



Au(x,t) + Vp(x,t) 


= t(x,t) 


in 2r, 


V-u(x,t) 


= 


in <2r, 


u(x,t) 


= 


on Ej, 


u(x,0) 


= n(x,T) 


on 


p(x,0) 


= p(x,T) 


on 


f(x,0) 


= f(x, T) 


on 



Here u^(^:, t) is a given target velocity, v > is a cost or regularization parameter, and | . | denotes 
the Euclidean norm in R d . We assume that u d (x, t) is time-harmonic, i.e.: 

Ink 

u d (x, t) = u d (x) e im with co = for some k e Z. 

Then there is a time-periodic solution to the original control problem of the form 

u(x,t) = u(x) e im , p(x,t) = p(x) e im , t(x,t) = f(x) e imt , 

where u(x), and f(x) solve the following time-independent optimal control problem: Min- 
imize 

\ [ \u(x)-u d (x)\ 2 dx+^- [ \f(x)\ 2 dx 

subject to 

id)Mu(x) - Au(x) +Vp(x) = f(x) in 

V-u(x) =0 in a, 
u(x) =0 on r. 

Using appropriate finite element spaces V/, of dimension n and Qh of dimension m for u and p, 
respectively, and the same finite element space \ h for f as well, we obtain the following discrete 
version: Minimize 

l(u-uJ*M(u-uJ + ^fMf 

subject to 

(WMu + Ku- D T p = Mf, 
Du = 0. 

Here the real matrices M, K, and D are the mass matrix, representing the L 2 -inner product in 
V/,, the discretized negative vector Laplacian, and the discretized divergence, respectively. As 
before, underlined quantities denote the coefficient vectors of finite element functions relative to 
a basis. 
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Completely similar to the discussion in the previous subsection we obtain the following re- 
duced optimality system, again after eliminating the control, here f. 



M 



K + icoM 
-D 






D r 





K 



ifflM 
D 



v 




v 



D r 









u 








p 









w 









r 








This system was discussed in If28l for the special case (0 = 0, which corresponds to the steady- 
state version of the control problem. The preconditioner constructed in [28J leads to convergence 
rates of the preconditioned MINRES method that do not deteriorate if the mesh size h and/or the 
cost parameter v approach 0. 

We will now apply the theoretical findings of the preceding sections and construct a precon- 
ditioner that will also work for general co. Convergence results will be derived that guarantee a 
bound for the number of iterations which is independent of h, v, and 00. If applied to the special 
case (O = 0, this bound is more accurate than the bound derived in ESI . 

A first and essential observation is that, by swapping the second and the third rows and 
columns, we obtain a saddle point matrix with a vanishing 2-by-2 block in the right lower part. 
This leads with a simple scaling to the system 



M v/v(K-ifl)M) 
v^(K + i(oM) -M 
-VvD 
vD 





vD r 






■VvD r 








u 








n/v — 









E 

















which is of the form (OQ) with 
A 



M 

y/v(K + ia)M) 



Vv(K-ifi)M)' 
-M 



-i2nx2n 



B 



-VV 



D 
D 



2mx2n 



and C = 0. 



(40) 



(41) 



Using the results of Sect. [3] and 0] we obtain the following properties of the bilinear forms a and 
b associated with A and B, respectively, in terms of the quantities a, \\a\\, j8, and \\b\\: 

Theorem 6. Let Jt be given by (0), and rt?71) . Then, for 3?, given by 



P 
R 



with P 



P 
P 



and R = v 



S 
S 



with P and S real and symmetric positive definite, the spectrum of 2? 1 j$ is real and symmetric 
around zero. Moreover, for 



P = M + v / v(K + (oM) and S = D [M + Vv (K+ Q>M)] D 



(42) 



22 



the following estimates hold 



a>^, NI<1, JS = 1, 



1. 



Proof. The symmetry of the spectrum around zero follows from Theorem[3] applied to the system 
in the original ordering of the rows and columns. Let 



H 



I (l-j)l' 

-i 



Then, by direct calculations, one shows that 

Re(Au,Hu) = \\u\\ P and = 3 ||u||p for all u G C 2 '\ 

Furthermore, observe that 



Therefore, 



Hu G ker5 for all u G ker£. 



• r |(A«,v)| . r |(Am,J/m)| 

a = inf sup „'\, ,', > inf 1 



O^GkerfiO^vekerS IMIpIMIp O^uekeiB \\ u \\p\\Hu\\p 

^ . Re (Am, //a) . Hull* 1 

> inf — = inf 



O^aekerfi IM|/ , ll# M l|/> O^MGkerfi || w||/>\/3 1| w||/> V3 



Moreover, since 



and 



(Au,u) = (Mu,u) 

+ 2 v / v(Re(Ku,w) - co Im(Mu,w)) - (Mw,w) with u 



|Im(Mu,w)| < |(Mu,w)| < ||u||m||w||m) 
|Re(Ku,w)| < |(Ku,w)| < ||u||k||w||k, 



it follows that 



\(Au,u}\ < ||u||M + 2 V / v(||u||K||w||K+£t»||u||M||w|| M ) + ||w| 



M 



< ||fi|lM+Vv(|N||+||w||| + fl) (||U|&+ ||w||y) + ||w||^ 

= (Pu,u) forallMGC 2 ", 

which implies that ||a|| < 1. 

Finally, since (BP- l B T q 1 q) = (Rq,q) for all q G C 2m , it directly follows from ([28]) and (|29l) 
that/3 = ||&|| =1. - - - - ^ 
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If a, || a ||, /3, and \\b\\ in (PT71) and (12TT) are replaced by the corresponding lower or upper 
bounds provided by Theorem [6J it immediately follows from (|24l) that the spectrum of 3P~ X 
is contained in the set 



--(l+V5),-/i 3 



U 



)U3,-(1 + V5 



where /I3 is the smallest positive root of the cubic equation 

These intervals read in 3 -digit accuracy 

[-1.618, -0.306] U [0.306, 1.618]. 
From the second part of Theorem \T\ we know a simple lower bound for /I3: 

1 



M3 > 



2^3 



0.289. 



Remark 2. The essential step for estimating the inf-sup constant a from below was the existence 
of the matrix H. Translated into a general Hilbert space setting and slightly more general, the 
essential requirements on such a linear operator H : V — > V are that 

1. there are some positive constants c\, C2 such that 

\a(u,Hu)\ > c\ ||w||y and ||.ffw||y < C2||w||y for all u G V, 

and 

2. ker5 is an invariant subspace of H. 

Then it follows analogously to the previous proof that a > c\jci. 



5.3 A numerical example 

We refer to ifTTl for numerical experiments for time-periodic parabolic optimal control problems. 
Here we present some numerical experiments for the time-periodic Stokes control problem on 
the unit square domain Q. = (0, 1) x (0, 1) cl 2 . Following Example 1 in lfT4l we choose the 
target velocity u d (x,y) = [(U (x,y),V(x, y)] T , given by 

U(x,y) = 10^-((p(x)<p(y)) and V(x,y) = -10^(<p(x)(p(y)) 

with 

<p{z)= (l-cos(0.8^z))(l-z) 2 . 
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The velocity Ud(x,y) is divergence free. Note that, contrary to the velocity tracking problem for 
time-periodic Stokes flow considered here, in [14] the velocity tracking problem was discussed 
for a time-dependent Navier-Stokes flow. 

The problem was discretized by the Taylor-Hood pair of finite element spaces consisting 
of continuous piecewise quadratic polynomials for the velocity u(x,y) and the force f(x,y) and 
continuous piecewise linear polynomials for the pressure p(x,y) on a triangulation of Q.. The 
initial mesh contains four triangles obtained by connecting the two diagonals. The final mesh 
was constructed by applying £ uniform refinement steps to the initial mesh, leading to a mesh 
size h = 2 . The total number of unknowns on the finest level £ = 4 is 18 056. 

Table Q] - [3] contain the numerical results produced by the preconditioned MINRES method 
with the block diagonal preconditioner as described in Theorem [6j The considered values 
for the mesh size h, the frequency CO, and the regularization parameter v are specified in the 
table caption and the first columns. The second columns show the minimal intervals [^3,^4] that 
enclose all positive eigenvalues of ^ . These intervals were computed with an extended 
version of the preconditioned MINRES method by the help Ritz values and harmonic Ritz val- 
ues, see lfl~6ll. The third columns show the (constant) interval [^3,^4] containing all positive 
eigenvalues as discussed right after Theorem [6] for comparison. The fourth columns contain the 
number k of MINRES iterations that are required for reducing the initial residual in the - 
norm by a factor of £ = 10 with initial guess xq = 0. This number k is compared with the 
theoretical bound k, which is shown in the last columns, based on the estimate 

IMI^^INI^ wtth g = — ^_ T , 

for the residual of the k-th iterate, see, e.g., [13], and the eigenvalue estimates for i^ -1 
discussed right after Theorem [6[ 



Table 1: V = 1, CD = 1 



h 




[jU 3 ,M4] 


k 


k 


1 


[0.627 , 1.595] 


[0.306, 1.618] 


6 


102 


0.5 


[0.620, 1.612] 


[0.306, 1.618] 


26 


102 


0.25 


[0.619, 1.616] 


[0.306, 1.618] 


28 


102 


0.125 


[0.618 , 1.618] 


[0.306, 1.618] 


28 


102 


0.0625 


[0.618 , 1.618] 


[0.306, 1.618] 


28 


102 



As expected from the results of Theorem [6[ the eigenvalues of ^ \M are bounded away 
from and 00 independent of h, CO, and v leading to a uniform bound for the number of MINRES- 
iterations. The theoretical bound 114 for the largest positive eigenvalues is rather close to the 
observed values, while the theoretical (uniform) bound /I3 for the smallest positive eigenvalues 
underestimates the observed values. Numerical experiments indicate that this is mainly due to the 
lower bound 1 / \^3, which underestimates a, and not due to the subsequent use of the estimates 
from Theorem [T] 
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Table 2: h = 0.0625, v = 1 



CO 




[M3,M4] 




k 





[0.618 , 1.618] 


[0.306, 1.618] 


18 


102 


1 


r / \ y 1 o 1 /; i oi 

L0.618 , 1.618] 


TA ""> A /l 1 /T1 oi 

[0.306, 1.618] 


O O 

28 


1 AO 

102 


1 A 2 

10 


L0.611 , 1.613] 


TA O A/1 1 /;i oi 

L0.306, 1.618] 


A O 

42 


1 AO 

102 


10 


rA C/1C 1 H A Al 

L0.565 , 1.614] 


TA O A /I 1 /T1 OI 

L0.306, 1.618] 


/I /I 

44 


1 AO 

102 


1 a8 

10 


rA /" 1 n 1 ^ i o i 

[0.618 , 1.618] 


rA o A z' 1 y 1 on 

[0.306, l .618] 


16 


1 A^ 

102 




Table 3: h = 


= 0.0625, CO = 1 






V 




[M3,/*4] 


Jfc 


k 


io- 8 


[0.566, 1.614] 


[0.306, 1.618] 


43 


102 


io- 4 


[0.611 , 1.613] 


[0.306, 1.618] 


42 


102 


10~ 2 


[0.619, 1.618] 


[0.306, 1.618] 


38 


102 


1 


[0.618 , 1.618] 


[0.306, 1.618] 


28 


102 


IO 8 


[0.618 , 1.618] 


[0.306, 1.618] 


28 


102 



6 Concluding Remarks 

The results from Sect. 2-4 apply to quite general classes of saddle point problems and provide 
sharp stability estimates as well as particular spectral properties. The subsequent discussion of 
two model problems from optimal control does not only demonstrate the applicability of the 
theoretical results but also shows the robust behavior of the associated preconditioned MINRES 
method with respect to the involved numerical and model parameters for the particularly chosen 
preconditioners. 

Concerning the implementation of the iterative methods one important issue has not been 
addressed so far. The use of a preconditioner & requires the evaluation of expressions of the 
form £P~ l y for some given vector v. In all discussed optimality systems for the discretized model 
problems this is a nontrivial task due to the potentially high number of involved unknowns. 
For example, the application of the preconditioner for the velocity tracking problem requires 
the evaluation of P ! v and S q with P and S, given by (l42b . for some vectors v and q. In a 
first step of approximation, is replaced by (l + y/v(0)M~ l + y/vcoK p l (Cahout-Charbard 
preconditioner, see 0), where M p and K p are the mass matrix and the discretized negative 
Laplacian in the finite element space Q/, for the pressure. With this replacement the application 
of the preconditioner involves only matrices which can be interpreted as discretized diffusion- 
reaction operators of second order. In a second step of approximation these matrices are replaced 
by efficient preconditioners (like multigrid preconditioners) which are well-established for this 
class of problems. Such a modified and efficiently realizable preconditioner leads to similar 
performance results as the original theoretical preconditioner according to the analysis presented, 
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e.g., in im @, ED, ED, ED. 



Appendix 

The detailed arguments for the inequality (|36l) . see the proof of Theorem[2l are as follows. 

The matrix A^ ax .P — A is positive semidefinite, see (1271) . Therefore, its Schur complement is 
positive semidefinite, too: 



^maxA-Aii-A 10 (A max / 5 o-Aoo) A 0l ^j u u u^ >0. 
Similar as in the proof of Theorem [2] one obtains 

\A W ^A maxJ P -Aoo) A ^j uuuxj > ^ A (AioA^Aoiwi, «i) . 



This implies 

(^max^l "All) Ml, Ml) - . A " (A 10 A^A 01 U h Ml) > 0, 

and, therefore, after dividing by (P\U\,u\), 

*L*-T2- 7A a _„ n>o (43) 

A max w 

for the Rayleigh quotients r\ and r 2 from the proof of Theorem^ 

Moreover, the matrix A — X^ m P is positive semidefinite. Then it follows analogously, how- 
ever only for A^ in < 0, that 

r 2-^L--—j^n>o. (44) 

a A min 

In the proof of Theorem [2] we need an upper bound for 

<Kn,r 2 ) = — ^— n-r 2 . 

It is easy to see that the maximum of the linear function 0(ri,r2) under the linear constraints 
(1431) and (1441) is attained at the intersection of the corresponding straight lines. By elementary 
calculations it follows that the value of (r\ , r 2 ) at the point of intersection of these lines is given 

by 

(^max ~l~ ^min ~ ~ ^max^min 

a- /I 

which completes the proof of (1361) . 
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